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This document contains my detailed calculation of the Generalised Few-cycle Envelope Approximation (GFEA) 
propagation equation reported and used in Phys. Rev. A67, 023813 (2003) [l( and its associated longer version at 
arXiv.org 0] ■ This GFEA propagation equation is intended to be applicable to optical pulses only a few cycles long, 
a regime where the standard Slowly Varying Envelope Approximation (SVEA) fails. 

The calculation is intended to be as complete as possible, but is still a "work in progress" , and so may, despite my 
qq ] best efforts, contain occaisional mistakes. It is an edited version of a longer document from which on-going work has 
f^i ■ been excised. Please contact me if you have any comments, corrections or queries. 

[*] I worked at this institution while doing the bulk of this calculation. My main project was with Prof. G.H.C. 
■ New on few-cycle optical pulses, and I was funded with money from the EPSRC. 

^ '. IMPORTANT NOTE: this calculation is of historical interest only: the current 
j-t , ! state of the art is summarized in larXiv:0707. 09821 Please note that the approach 
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Nonlinearity: B(rj_,£, t; A) approximatio: 



I. INTRODUCTION 



This report is a calculation intended to generate an evolution equation for an envelope approximation description 
of pulse propagation in the few-cycle regime. I try to fully investigate each step and approximation, include all the 
algebra, and discuss any subtleties that may arise. 

The calculation parallels and extends that in -kit T. Brabec, F. Krausz, "Nonlinear optical pulse propagation in 
the single-cycle regime" Q. I also (amongst other things) make some comments about their slowly-evolving- wave- 
approximation (SEWA). 

The calculation includes all the steps taken from the Brabec and Krausz Q starting point of (their eqn.(l)) 



(dl + Vi) E(r, t) - 1# f dt'e(t t')E(r, t') = %d 2 t P(r, t), (1.1) 

c J — oo c 

through to their basic nonlinear envelope equation (NEE) (their eqn.(6)) 



dtA = —£A + iDA + — — 1 + — <9 T VIA + i—^- 1 + — <9 T B. (1.2) 
2 2/3 V w ) n o \ w o / 

This report also derives the NEE equation found by M.A. Porras, "Propagation of single-cycle pulsed light beams 
in dispersive media" [3] which extends the Brabec and Krausz Q theory so as to treat diffraction and self- focussing 
better, but which neglects the nonlinearity, giving their SEEA equation: 



= -Y A + M + Wo{ 1 + ^) ViA (L3) 

There seem to be hints of the Brabec and KrauszQ result in the early paper by + J-A. Fleck, "Ultra-short pulse 
generation by Q-sqitched lasers" [5|. Fleck[5[ has an envelope eqn.(2.3a,b) of (not using eqn.(2.2)) 



c 

-d t £- - d z £- + 2ik = - (1 + iS) P' + c.c, (1.5) 



±d t £ + + d z £ + = - (1 + iS) P' + c.c, (1.4) 
c 

V. 



where the 1 + iS results from a 1 + d±/u> term acting on the polarization, which has a linearly decaying memory 
with characteristic time T2 (and i5 = (0JT2) 1 ); dispersion and diffraction are not included. It is seems plausible that 
we can get from eqn (ll.2p to (|1.5[) - we scale the variables (the dt will disappear with a co- moving frame), neglect 
diffraction, neglect dispersion, and introduce the same polarization model. Etc. This will be clearer as we follow the 
full calculation in the next sections. 

An early review of attempts to describe ultra-short pulse propagation was given by Lamb Q . 
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II. ENVELOPES AND CARRIERS 

The substitution used in the "envelope approximation" is just a splitting of a general waveform (of e.g. the electric 
field amplitude) into two parts, an "envelope" , and a "carrier" ; where the carrier part is intended to carry almost all 
of the oscilliatory part of the waveform. It is usual to describe a plane-polarised wave, with the B field perpendicular 
to, in phase with, and proportional to (in amplitude) the E field. It is also possible to do such a calculation with 
circularly polarized field variables 0] (using E c — E x + iE y rather than just E x ). 

For a forwardly propagating carrier exp [i {(3qz — u>ot + ■f/'o)] an d an envelope A(f±, z, t), the substitution is 

E(f,t) = A(r ± ,z,ty iPoZ ~ Uot+M + A^r ± ,z,t)e-^ oZ -" ot+M (2.1) 

= A(r ± ,z,t)e lS ~ + A*{r ± ,z,t)e~ lS ~ , (2.2) 

^ = (ftzTwot + W- (2.3) 

Note that the carrier is forward propagating because of the chosen signs on the wavevector (/3q) and frequency (loq) 
parts: both terms in eqn. (|2.2p are forward propagating (c.f. a wave described by e.g. f(x — vt)). Also note that the 
Poynting vector for the carrier field also has a direction, given by E x B - and here B is determined by E (c.f. FleckQ 's 
approach). Using a forward carrier means that any backward propagating components that happen to be in the E 
waveform need to be contained in the envelope - unless extra backward carrier terms are added to the substitution 
above. However, since in many useful situations an initially forwardly propagating wave does not develope a significant 
backward propagating component, we can use approximations to work in a regime where backward contributions are 
negligible, rather than complicate our representation of the E waveform: a forward+backward carrier equation might 
look like this: 

E(r,t) = A{r u z,t)e l{PoZ - u ° t+ ^ + A*(r_ L ,z,i)e- 4(/3oZ -"°* +v '+ ) 

+B{r u z,ty {f3aZ+UJat+ ' l) - ) +B*{r^z,t)e- l ^ aZ+ulot+,l '~ ) . (2.4) 

Typically we then try to factor the carrier part out of our equations of motion for the waveform, and simplify the 
envelope equation of motion by making approximations based on assumptions of (e.g.) smoothness of the envelope, the 
small contributions from backward propagating terms, and so on. This then leaves us with a (hopefully) manageable 
equation for just the envelope function. Note that our choice of carrier frequency is only constrained by the need to 
keep the approximations manageable. 

One point with the use of the envelope function substitution is that any phase- like properties of the waveform become 
obscured. This is because they can be contained in either (or both) the carrier and the envelope. An alternative 
choice in the phase of the carrier (e.g. replacing ipo with some tp' ^ ipo) will in fact mean that the envelope function 
is different. Plotted on a graph, this can seem to have a large effect, particularly for a few-cycle pulse, where the peak 
of the waveform moves about noticeably if the phase of the carrier is changed in a envelope-carrier pair - see fig. 1 
and comments in Brabec and KrauszQ for some analysis. 

Note however that the envelope (and carrier) are complex - and a phase shift of exp [i (ij)' — ipo)] m the carrier 
can be exactly matched by a fixed shift of exp [—i (ip' — ipo)] in the envelope. This means there is no calculational or 
simulation problem associated with shifting the carrier phase, as long as the envelope is also adjusted. What matters 
is that the real E field waveform resulting from a envelope-carrier pair is correct according to the boundary/initial 
conditions. Confusion can only occur from the point of view of an envelope-only picture of the pulse shape. 

In the usual case of a pulse that is long compared to its natural carrier frequency, the envelope should be smooth and 
so the relationship between the underlying carrier phase and the electric field EQ is often ignored, or, if not ignored, 
then can be regarded as specified by the many oscillations in E at the carrier frequency. This tend to lead one towards 
an "envelope only" view of the pulse, which can later cause confusion for few-cycle pulses where a different choice of 
carrier phase leads to a different looking envelope on a E vs z (or t) graph, where only the real part is plotted. For 
some kind of complex plot, we would clearly see that these different looking envelopes arc complex rotations of each 
other, and their "differences" are merely an artifact (but see the following subsection). 
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FIG. 1: F-envcar: Diagram showing how a field Ei is decomposed into an envelope A4 and a carrier H, the envelope A is 
propagated to its final state Af, then the final field Ef is reconstructed. It also shows that an alternate field E[ can have the 
same initial envelope when a different carrier H' is used, and that the same propagation can be used to extract E^s final state 
E'f 

A. A Phase Function 

Some authors (e.g. @, use a separate phase function 2p(r,t) in order to ensure their envelope function remains 
real. In this case the definition in ean. (|2.3|) looks like 

E(r,t) = Aflfrl.z,^ 1 ^ 2 ^ '^^'*' +A R (f ± ,z,t)e- t{f3oZ ~ UJot+ ^ t)) . (2.5) 

Whilst this might look like a good idea, it complicates any kind of propagation equation for the pulse envelope that 
we might derive: it would now contain additional derivative terms (of the phase function) and also we would need a 
propagation equation for the phase function itself. Further, this phase function is ambiguous (recall tp + 9 is the same 
angle as 9), would be undefined (or any value) when An is zero, and numerically difficult to handle when Ar is small. 

I do not use this sort of phase function in this document. 



B. A Single Envelope Represents a Set of Pulses 

For a given set of initial conditions (usually just the electric field profile of the input pulse(s), we might pick any 
value of carrier phase we liked. Each different value of carrier phase would result in a different pulse envelope - but 
each of the resulting combinations of carrier phase and envelope will specify the same initial conditions. 

So, starting with a fixed field, you can use a variety of carrier choices, and end up with a variety of envelopes. After 
solving for the propagation of the chosen "initial state" envelope, we get a "final state" envelope. From this we can 
reconstruct a unique final state E field. This parallels the left hans side of Fig. [T] 

But note since a given envelope may be turned back into an electric field by applying any values for carrier phase, 
this final state envelope can be used to generate a range of final state electric fields. Each of these fields corresponds 
to the initial condition specified by the initial state envelope and that same choice of carrier phase. 

So one envelope simulation provides a range of Ei — > Et solutions, as indicated diagramatically on fig. [T] 

Note that for polarization terms with their own dynamics (e.g. a two level atom, see my report 
Two level atoms and the few cycle regime |9| ) , the choice of carrier phase alters not only the pulse envelopes but 
also the representation of the initial polarization state. 
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FIG. 2: Diagram showing traditional "carrier phase" estimation method. The electric field is plotted a sawtooth because it 
was quick to draw - a modulated sinusoid would be a more realistic profile. 



C. The so-called "carrier phase", i.e. the pulse phase 



Many authors publishing on short pulses refer to "the carrier phase" . By this they seem to mean something derived 
by comparing an inferred field envelope to the peaks and troughs of the actual (oscillating) electric field. For example, 
they might take the distance between the peak of the envelope as their zero (i.e. the point of reference), and the peak 
of the nearest electric field oscillation as giving the field phase (see fig . A good example of this is the paper by 
Chelkowski and Bandrauk [lfj . To work, this "maximum amplitude" procedure assumes a number of things: 

(al) The envelope has a single peak: Of course this is not always true. Further, even if it does happen in some 
particular case to be true, that peak is not always well localised - such as in flat-topped pulses. 

(a2) The pulse does not have a complex phase structure: Of course it is straightforward to generate real val- 
ued envelopes from an oscillating electric field profile, as the electric field is real valued. However, if (e.g.) the 
pulse is chirped, the electric field oscillations will no longer have a fixed frequency, so it is not possible to use 
and envelope-peak to field-peak distance to guess a phase without additional assumptions. 

I believe it is unhelpful to talk of "carrier phase" in the contexts the term is usually applied. From a mathematical 
point of view, there are two phases: the carrier phase, and the envelope phase; further, the envelope phase may well 
have a complicated structure that obscures or overrides the role of the carrier phase. The carrier phase should be 
fixed according to some spatio-temporal reference point, then further discussions along the lines of "the phase of the 
pulse" should refer to the envelope phase, and so the phase structure of the pulse will not so easily be ignored or 
sidestepped. However, in the event that a single phase parameter is appropriate or desirable, a clearly defined method 
should be used to extract it, avoiding unreliable and subjective arguments about which value of the phase profile of 
the envelope is "the phase". 



III. EXTENDING BRABEC AND KRAUSZ: THE POST-TRANSFORM ENVELOPE 

Brabec and KrauszQ consider the case of small transverse inhomogeneities of the polarization, and so start with 
the three dimensional wave equation 

(dl + Vi) E(f, t) - ±d? f dt'e(t - t')E(f, = % (d t © 9 -d 7 ) d t P nl (r, t). (3.1) 
c J—oo c \ n J 

Here (as in Brabec and KrauszQ) is the transverse Laplace operator, d a is used as a shorthand notation for 
d/da, e(t) — (27r) _1 f_ dtae(u})e luJt , e(uo) = 1 + 4irx{u), and x( w ) is f ne linear electric susceptibility. The electric 
field E propagates along the z direction. Both E and the nonlinear polarization P n i are polarized parallel to the x 
axis. 



G 



The @ is a new term include to allow the calculation to apply to the G Fleck|5j] variables. For normal cases 
set g = and forget it. Here I use an alternative "±" sign: © - this is to distinguish this variable sign from a later 
(and independent) variable sign ± cause by the choice of carrier direction. 



A. The linear electric susceptibility 

Now we need to treat the effects of the linear electric susceptibility. We start by fourier transforming the equation 
using exp (— iu>t); but neglecting to keep track of the normalisation, since this will take care of itself when we transform 
back. Using the correspondance <9 t <-» — ioj, the first LH term is simple, and transforms to: 

(d 2 z +Vi)E(r ± ,z,u). (3.2) 
The RH term is also simple, and transforms to: 

^(- w @3Cq\ (- luJ )P nl (r^z,u:). (3.3) 
The second LH term (with the dt' integral) is more complicated: 



+ oo i pt 

dt e" lult -^d 2 / dt' e(t - t')E(f x , z, t') (3.4) 



(-lujf 



dt e~ wt / dt' e(t - t')E{f ± ,z, t'). (3.5) 



If the upper limit of the dt' integral was t and not oo, the dt' integral part would be the normal convolution integral; 
hence we could convert it into the product of the fourier transforms of its constituents. This could be justified by 
saying e must be causal, and so is for any t' > t, hence the limits of the integral can be extended to +oo. 

Extending the upper limit of the dt' integral to oo gives 



-lUjt 



u 2 „, 



— oo 



— / dte~ %wt / dt' e(t-t')E(r ± ,z,t') (3.6) 

C J — J co 

dt e- wt e(t)\ dt e- wt E{r ± , Zl t)^ (3.7) 

+^f»E(r 1 ,z, W ). (3.8) 
The resulting transformed version of eqn. (|3.1|) is 



(d 2 z +Vl) E(r ± ,z,u>) + ^-e{u>)E(r ± ,z,uj) = (-vj © 9 -^d z ) (-tw) P n i(?±,z,Lj). (3.9) 

Now I might want to expand i(u>) in powers of uj, but to make things easier I'll replace it with the fc 2 (w) and expand 
k about uj e instead. Using e(ui) — c 2 k{uj) 2 / lo 2 (as do Brabec and Krausz[3|) and then 



k{uo) = ~\ ; 7 " = 9 - fc(w)l ^ + /3n,a n eR; a n = Im( 7 „), /3 n = Re( 7 „). (3.10) 

Note that Brabec and Krausz[3l have a n /2 in their equations where my definitions will give a n . This is because 
the definition Brabec and Krausz[3j give for a n (below their eqn. (BK3)) is not the one they actually use. The one 
they use is consistent with a corresponding to the decay in the intensity, not the field; and in fact PorrasQ alters 
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his definition of a n from that stated by Brabec and Krausz[3j (and me) in order to have terms like a n /2 appear. 
Porras 1 . despite his different definition, is consistent with my calculations. The a n that Brabec and Krausz[3j use 
is the same as that defined by PorrasQ; and both Brabec and KrauszQ and Porras [H use uj e = loq. 
Using this expansion, the equation becomes 

2 2 ^ 

(ft 2 +Vi) E(r ± ,z,u;) + ^^k(Lj) 2 E(r ±7 z,u;) = (-iu> © 9 — ft) {-icj) P(r ± ,z, w) (3.11) 

c ui c . V n I 

This can then be transformed back into the time domain (NB: ft <-> — id, id t <-> u>; (w — — » (ift + i 2 lo c ) = 
i(d t + lU) e )) 



(d 2 + V 2 ± )E(f ± ,z,t) + 



B. The envelope and carrier 

Now I split the field up into an envelope part and a forwardly propagating carrier-wave part using the substitution 

E(f,t) = A(r L ,z,t)e j(/3oZTa,ot+,/ ' o) + A*(f_ L ,z,i)e^(' 3oZ=F " ot+,/ ' ) (3.14) 
= A(f ± ,z,t)e lST + A*(f ± ,z,t)e- lST , (3.15) 

and similarly for P„;(r, t) — B(r±, z,t; A)e 1 ^ + B* (fj_, z,t; A)e~ l ^ f . The symbol S T is introduced purely as a 
convenient shorthand notation for the terms in the exponential; and the minus sign (i.e. S~) refers to a forwardly prop- 
agating carrier, and plus sign (i.e. S + ) a backwardly propagating carrier. With these envelope-carrier substitutions, 
the equation of motion becomes 



E 

n=0 



7„ (w - uj e y 



E( 



r^,z,LU) = 



—iuj fcB — a. 



(-lUj)Pnl(fx,Z,Uj). (3.12) 



E 

_n=0 



«"7n (ft + lUeY 



E(r ±7 z,t) = %(d t ®^d z \d t P n i{r±,z,t). (3.13) 



([lA, + d z } 2 + Vl) A(f x , z, t) + e- HT ([-fa + d z } 2 + Vl) z, i) 



E-T^ (T^o + ft + «w £ ) n 
77.' 



n=0 

' C 2 
^2 



A(f ± ,^,t) 



(±«^o + ft ~ JW e )™ 

77.' 



A*(r ± ,z,i) 



(T^o + ft) ©^(i/3 + ft) (^o + 9 t )S(r ± ,z,t;A) 

n J 

(± lWo + ft) © — (-»A) + ft)l (±K4> + ft) S*(rL, z, t; A), 

2)1 



(3.16) 



(NB: i n (tiuo + ft + -» (tw e T + ft)" - * 2 " (w e T^- *ft)" -> (-wo) B (T T - *d t /^ ) n ) 



([»& + ft] 2 + Vl) A(r x , «, t) + e- HT ([-»#, + ft] 2 + Vi) A* (r ± , z, t) 



oo / s n 



.ra=0 

^4tt 

' c 2 

72" 



w 



A(rj_,z, t) + e 



E 

.71 = 



7„ (-w ) n 



(T^o + ft) © — (t/% + ft)l (T^o + ft) 2 B(f ± , z, t; A) 
n J 



[±iu) + ft) 



0C 



H/3 + ft) (±KJo + &)B*(»*_L,«,t;A), 



A*(ri,z,t) 



(3.17) 



Here = (uj e =F ^o) /^Oi which is a quantity which usually would be set to zero - but retaining it allows me to 
expand the dispersion around a frequency other than ojq. Note the usage of T T and w € is clumsy, because we need to 
alter its sign under complex conjugation; carrier direction reversal is taken care of with the ± notation. 
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I now split eqn (|3.17[) into two separate equations, the first "A equation" containing the terms like e 1 " , and the 
second " A* equation" containing the terms like e _l ~ T . These two equations are simply the complex conjugates of one 
another, and so writing down only the first one is sufficient: 



e lST (jtfo + d z f + Vi) A(f ± , z, t) + e lS * 



30 



E 

,n=0 



7™ (~^o) r 



- —d t 

iv 



A(f x ,z,t) 



(=Rw + dt) 



(if3 + d z ) \ (Tiuo + dt)B(rx,z,t;A). 



(3.18) 
(3.19) 



This is simplified with number of minor steps: dividing by the e 1 " factors (which are conveniently never zero), 
extracting factors of =F(t>o from the RHS ( => (^ftuo + d t ) = (=R w o T v^odt / (tm^o)) = T l ^o(l + <9 t /(=Fiwo) = T iW o(l T 
dt/(uiio) = ^fitOoO- i*^t/wo), then preparing to use 1 = c/3o/nuo, leaves 



(b(3 + d z } 2 + V'i) A(f x ,z,t) + 

47TWq 



E 



7n (~^o) ? 



1±— 5 t 

CJ 



,n=0 

nw V Po 



-i 2 



A(f±,z,t) 



l±—d t B(rx,z,t;A). 

U 1 



(3.20) 
(3.21) 



This appears to differ slightly from Bracbec-Krausz eqn. (2) in that it has the opposite sign on the RHS - however, 
agreement is recovered later in eqn. (|3.48p . 

If Aq is a solution of the A equation, then its conjugate Aq is a solution of the conjugate A* equation. This means 
that solving one solves the other, and a total waveform can then be easily reconstructed using eqn (|3.15p . There are 
no approximations made in doing this, but it may be that there are complicated (or subtle) cases where solutions of 
the full equation are not expressible in terms of solutions of the two separate ones. 

Note that for the forward propagating carrier (S - , the upper sign choice in S^), there are no explicitly backward 
propagating terms, despite the fact that we have not excluded them in any way - this is because they do not 
arise spontaneously, but need to be created. Such an effect could occur in the case of multi-field systems or exotic 
polarization behaviour, where there may well be spatially oscillating terms from the nonlinear polarization term B 
( e.g. exp {±.iksz)) which could force a (possibly backward propagating) oscillation onto A. This would very likely 
violate some approximation we will want to make later, e.g. a "smooth" or slowly varying envelope function A. We 
can only neglect backward propagating components if (a) there were none to start with, and (b) by verifying (or 
assuming) that the nonlinear polarization has convenient properties - although we could extend eqn (|3.15p to include 
backward carrier terms as already discussed (see section [TTJ) . 

At the equivalent point to eqn. (|3.2ip in the Bracbec-Krausz paper, they already claim to have neglected backward 
propagating waves (after their eqn (1)): "the neglect of backward propagating waves is consistent with the approx- 
imations that will be made in the following derivation of the envelope equation and will be commented on later" 
- their comment being that "excessive" change in the envelope can lead to backwardly propagating components to 
the envelope (see e.g. Shen [HI]). ■"•> however, leave any approximations relating to the neglect of backward terms 
to later on - it is still the case that (in principle) that the envelope function might contain backwardly propagating 
components. 

This differs from Fleck H in that his E ± and £ ± (G+ and Q+ in my notation) are constructed as explicitly forward 
and backward propagating, and hence should not really be compared directly to my E or A. My E is in fact Fleck[5|'s 
E + + E~ (G + + G~), and somewhere in the approximations used to get to my starting point eqn. (|3.1|) the magnetic 
field parts (retained by Fleck Q) have been assumed to be solely dependent on (deriveable from) the electric field. See 
my derivation "A second-order wave equation using Fleck field variables "fl2l. fX3(] for more information. 



C. Scaled co-moving variables 



I now change into a scaled co-moving reference frame, but one slightly different to that of Bracbec-Krausz - the 
difference being that I scale space by 0o and the time by ujq as well as shifting the origin. Bracbec-Krausz and Porros 
use £ = z and r = t =F 0iz, but I instead put 



T = LJQ (t T P[z) . 



(3.22) 
(3.23) 



Here I use (3' and f3[ rather than just po and j3\ because it may not always be convenient to use the natural scaled 
co-moving reference frame; which might well be the case for solving multi-mode problems. For an A* equation, the 



signs of both z and t should be reversed so that = — [3' z and r_ = 
frame (for S~ , upper sign) is backwardly propagating, since although 
derivatives for the A equation transform like 



-u> {t =F P[z). Do not think that this A* 
—z, time has also reversed t_ ~ —t. The 



d t = 

d z = 



d_ 
~dt 
d_ 

dz 



dr d ^ g?£ d d _ ^ 

dt dr dt d£ dr 



d£ d dr d 
dz dt; dz dr 



a d qi d 



(3.24) 
(3.25) 



(for A*, d t — —oj d T _ and d z — —P' d^_ ± uJo^dr ). The scaled co-moving A equation is then (with q — (3' /po). 
Note that now a condition that £ 3> 1 refers to lengths 3> 1/po (i.e. "long"), and t>1 refers to times 3> l/w (i.e. 
"slow"). Similarly, £ -C 1 refers to lengths -C l//3o (i-e. "short"), and t<1 refers to times Cl/wo (i.e. "fast"). 



(3.26) 



n=0 



n) 



+ - 



(1 ± ld T ) 



gcf3 



(1 — iqd^ — ia'd T ) 



+ ^a2 [(^ © l 9 a 'd T ) © 5 (1 — *<Z^)] (1 ± i<9 T ) B(f±, £, t; A) ... now expand, prepare for £>' .. 

c p 



/?o 2 



(1 ± i<9 T ) B(r±, £, r; A) ... now divide through by pg 
2 ' 

' A(fj.,e,T) 



(3.27) 



oo / % n 



,n=0 



Po 



A) 



po ia wo/?i / t t ~ n . w ~ 



47T 

+ — [1 ± z9 T © t 5t 7'3 T © g (1 - i«0 £ )] (1 ± id T ) B(? ± ,t, r; A) 



Po 



"oPi 



1 + 2i [qds T ^f±d T ) + [qd ( T ~^d T + -y V± 



Po 



PI 



+ — [1 ± id T © igo-'d T © g (1 - «j0 £ )] (1 ± id T ) B(f ± ,£, r; A) 



Po Po 
... and now divide by 2i... 



(3.28) 
(3.29) 
(3.30) 
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u_oP[ 
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2i 
2~i 



'^oP 1(T T_ ldr) + :- _ d y 
Po 
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I 2 


' UoP\ 


+ 27 


1 Po 



vi 



- (T* - id T ) 



2ir 

+ — [1 ± id T © z. 9( r'a T © .9 (1 - «f0 £ )] (1 ± id T ) B(? ± ,t, r; A) 



mf 



2i 



vi 



ZCT (T T — ld T ) 



Po 



Po 



1 2 



ZD' 



ZD' 



^(r±,£,r) 



(3.31) 



za (TT - td T ) + 

Po 



n 2 



^(r±,£,r) 
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2tt 

2 [1 ± id T © iga d T 

in 



,(l-iqd $ )}(l±id T )B(r ± ,Z,T;A). 



(3.32) 



Here I have introduced the dimensionless a = woA/A = (wo/A)/(1/Al) = v fl v gi <*' — w oAi/Aj an d used the 
fact that the refractive index at ujq is uq = cA/^o- I also define a dispersion term D in a similar way to Brabec and 
Krausz [Ij], but instead use a scaled (dimensionless) version D' in following equations: 



D' = ^D 

A 



wo 

~A 



iai (T T - id T ) + J2 



Tn (-^o) 



(T T - id r ) n 



n=2 



(3.33) 



with 7^ = 7„ for n > 2; otherwise 7q = 0; 7^ = tai: the parameters «o, Aj A are handled separately from D 
because of their important role. 



D. Aside: Brabec and Krausz approximation criteria 

Brabec and KrauszQ introduce some criteria designed to motivate approximations to their equations (BK5a,b,c). 
Since I use the same variable names, but differently scaled, I here write the Brabec and KrauszQ criteria in their 
form on the LHS, and indicate with an arrow my form on the RHS. 



(5a) 
(56) 

(5c) 



« A |^| 
\d T A\ « lo q \A\ 



\d(A\ « A, 
\d T A\ < A, 



A) 



(3.34) 
(3.35) 

(3.36) 



Note that g — here. The motivation for the first two (5a, b) are obvious from ecm. (|3.31[) . in my scaled co-moving 
frame - they allow me to say certain quantities are small, and hence I could choose to neglect them. In contrast, 
it is not clear how the third condition will make my equations simpler beyond removing a single prefactor (since 
|woA/AI ~ 1); but |woA/A -C 1 might seem better still, if perhaps not as physically relevant. However, Brabec 
and Krausz [3] collect their terms together differently, and indeed the situation becomes clearer after I rearrange the 
eqations. 



E. The Generalised Few-Cycle Envelope equation 

I keep (=pT — id T ) terms intact because for the usual case where ui e — luq is chosen, they will simplify to — id T . Note 
also I have divided the equation through by Aj rather than the single A of Brabec and Krausz [3j. Still retaining all 
terms in the equation, I have 



= 



(qdz =F a'd T ) + i =F a'd T f + 



(T* - id T ) 



1 

+ 2 



(a) 



x A 



2tt 



[1 ± id T © iga'd T © g (1 - iqd t )\ (1 ± id T ) B(r x ,£, r; A) 



(qd e T o'd T ) + - (qd e =F & ' d T f + la (T* - td T ) + ( 2° - t& 
11 v po 



-a 2 (T* - id T y + 2ia (T=F - id T ) ( g - it)'\ + I ^ - 4j D' 



A) 



-tD' 



[1 ± id T © i .ga'a T © (1 - zoc) ? )] (1 ± id T ) B(f±,(, r; A) 



(3.37) 



VlU(ri,e,r) 



(3.38) 
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(a') = Uqdt T o>dr) + j t (qdz T °'d T f + la (T T - id T ) - '-a 2 (T* - % d T f 



*(T*-ia r )(£-ifl') + (£-*/y 



Po 



a 



0o 



2 \0o 



1 



2i/3$ 



2tt 



[1 ± id T © iga'd T © g (1 - 19^)] (1 ± *d r ) B(?x,£, r; A) 



(b) = \qd£T o'd T + jq 2 d 2 =F ^2^a'9 r + ^cr ,2 d 2 + iaT^ + ad T - ^- (T* - *d T ) 2 - cr (T T - z5 r ) 



2/ 



2? 



A) 



(3.39) 



-tD 7 



f£ - l b'\+ L \ ?L - 1& 
Po 



2 \0o 



ViU(fx,e,r) 



2tt 



in? 



[1 ± z9 r © iga'd T © g (1 - ^)] (1 ± id T ) B(r x ,£, r; A) 



(3.40) 



(c) 



|(1 ± za'a T ) ^ + z(T + 1 (q 2 2 + a' 2 d 2 ) + (a T a') d T - ^ - id T ) 2 - a (T=F - z5 r ) - 



2 V# 



viU(fx,e,-r) 



2tt 



[1 ± id T © zga'd T © g (1 - z^)] (1 ± id T ) B(f x ,t, t\ A) 



(3.41) 



(d) 



(1 ± za'cV) g0 £ + zcr + 1 (g 2 df + cr' 2 c> 2 ) + {a T </) d r - ^ (T T - z5 r ) 2 - a (T T - z5 r ) ^ - iU 



~2~ 2 



Po 



2t 

2 



2i/3( 



;ViU(rl,e,T) 



2tt 



[1 ± zd r © iga'd T © g (1 - (1 ± id T ) B{r x ,£, r; A). 



(3.42) 



This equation (|3.42|) has a nice 1 4- i(a>o//3o — tb') term which may be useful in a generic "small dispersion" case. 

However, we want to reach the more typical Brabec and Krausz|3f (or Porras|4|) form; and so instead continue with 
the algebra from eqn. (|3.41|l . Like Porras[4|, I retain the natural choice of 1 + iad T multiplier for the term. This 
differs from the 1 + id T used by Brabec and KrauszQ, and is why PorrasQ claims that the Brabec and KrauszQ 
SEWA equation has the space-time focussing "slightly falsified" , mentioning RothenburgpjJ ~ it is because Brabec 
and Krausz 3] had already introduced a a = 1 approximation, and as such the "missing" cr is an approximation, not 
anything "slightly falsified" . 

Now ruthlessly shift all the T=F, 0(<9f) and 0(<9 2 ) terms to the RHS... 



(1 ± ia'd T ) qd^ + (cr T cr') d T + (1 + iad T ) 



00 



l D')\A(r ±7 Cr) + ^V 2 ± A(r ±: tT) 



^ [1 ± id T © iga'dr © g (1 - tqfy)] (1 ± id T ) B(f x ,t, r; A) 



rT T - aT T 



(To 

00 



ib' 



2 V T > 2 \p 



it)' 



2i 



(q 2 d 2 + cr' 2 d 2 ) 



A(rl,£,r).(3.43) 



Presumably Brabec and Krausz Q used their eqn. (BK4) to motivate their (BK5c) cr ~ 1 to get a (1 + id T ) term 
- because they wanted to cancel the same factor from the nonlinear term multiplying B(r x , £, r; A). This can be 
contrasted with PorrasQ, who was not interested in the nonlinear case but wanted instead to cancel terms in the 
"V^" diffraction (or self- focusing) term. In fact using Brabec and Krausz Q eqn. (BK4) as a midpoint on the way 
to Brabec and KrauszQ eqn. (BK6) is not the best path, even for Brabec and KrauszQ - their early selection of a 
(1 + id T ) prefactor for is unnecessary. I show this here, but first hide all the T dependent terms in Ty (as they are 
usually zero anyway), and the rest of the RHS terms in Trhs because they remain unchanged. 

After some sign changes caused by moving factors of i from the denominator to the numerator, and by moving one 
inside a set of () brackets, we get 



12 



= (l±i<j'd T )qd(:A{r ± ,S,,T) + (crT<T')d T A(r ± ,£,T)-(l + iad T ) 



Oo 



+ iD')A(r ± ,tr) 



2nr 



\7 2 ± A(r ±7 t t) - — [1 ± id T © iga'd T © 5 (1 - iqd^)} (1 ± id T ) B(r ± ,£, r; A) - Ty - T RHS , (3.44) 



Trhs — 2 



-, 2 af + (^-<,' 2 )a; + 1^->d' 

1 Po 



(3.45) 
(3.46) 



Hence 



qd 5 A(r ± ,£,T) 



-d T A(ri_,£,T) + 



1 ± iC r'<9 r 



TZ — 77T [~lT + lD ) A ( r ±^> T ) 
1 ± ut'Ot \ po 



201 (1 ± icr'<9 T ) 



+ 



2m (l + zd T ) 
rag (l±icr'5 T ) 



B(r x ,t,r;A) + 



Ty + Tflgs 
l±z<j'd T ' 



ViA(r x ,£,r) 



(3.47) 



Now to avoid overly complex equations I set a 1 = ±<t (i.e. (3[ = p\), and a carrier appropriate to the group velocity 
by choosing the upper sign), forcing the r scaling for the field to match the material rather than (e.g.) another field. 
Note that in the case of multiple field components with different group velocities, it may be necessary to have f3[ ^ j3\ 
in order to keep all of the co-moving frames aligned. 

One final step now gives us a Generalised Few-Cycle Envelope equation for a propagating pulse in a nonlinear 
medium. It is in the style of Brabec and KrauszQ, but unlike that of Brabec and KrauszQ (and of Porras[J]), it has 
no approximations beyond that from the starting point, the permittivity convolution, and the separation of A and 
A* propagation equations. Remembering to expand the Trhs and TV terms as necessary, with g = 0, we have the 
Generalised Few-Cycle Envelope Approximation (GFEA) equation: 



qdzA(f±,£,T) 



Q'O 



^+iD') A(r ± ,^r) 
Po 



2/3g (1 + iad T ) 



2m (1 + id T ) 
n\ (1 + iadr) 



Ty + Trhs 
1 + lad-j 



(3.48) 



NOTE: If g ^ 0, replace one of the (1 + id T ) numerator terms multiplying B in eqn. (|3.48p with 
[1 ± id T © iga'd T © g (1 - iqd$)]. 



F. The nonlinear "few cycle" term 

We might prefer to handle the nonlinear term by basing it on (1 + id T ), for example if we intended to neglect the 
diffraction term entirely (e.g. as in Brabec and KrauszQ). By starting at eqn (13.48[) we see the possibility of a nice 
expansion in 5 = 1 — a by simply replacing the entire B prefactor term in eqn (13.48j) with: 



{i + id T y 



B(f ± ,tT;A) 



(1 + ld T ) [(1 + l(ld T ) + I (1 - a) d T ] 



(1 + ld T ) 



{l + iad T ) 
i (1 — a) d T 



B(f±,t,T] A). 



(l + iad T ) 

We could instead rearrange the B term by aiming at the 1 + iod T form of Porras 



(3.49) 
(3.50) 
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(i + id T y 

(1+M7<9 r ) 



B(f ± ,Z, T ;A) = 



[1 + iad T + i (1 - a) d T f 
(l + i<Td T ) 

(l + iad T ) 2 + 2{l+%ad T )i(l-a)d T - (1 - a) 2 d, 



(1 + iad T ) + 



(l + iad T ) 

2i (1 - a) d T + (-2cr + 2a 2 - 1 + 2a - a 2 ) d 2 



(1 + iad T ) 



l + (l-a 



{l + iad T ) 

2id T + (!+£) d 2 
(1 + iad T f 



(3.51) 

B(rj_,£,r;A) (3.52) 
B(fj.,e,r;A) (3.53) 
(3.54) 



It would be better, of course, to expand the few cycle term to a fixed order explicitly - there will be many expansions 
like the above, that when truncated to first order, are correct to within terms of second order; but which differ from 
each other by amounts that are also of second order. So: 



(l + tad T ) 



(l + 2id T -dl) x (l + iadrY 1 



(1 + 2id T - d 2 
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iad T 
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14 
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a)d T 



la 

iad T - 2i 2 ad 2 + iadl 



3 3l 



3! 

+ 0(0?)] 
a 2 d 2 - 2io 2 dl 



• o(a*) 



iu 6 d: 



O(0 T 4 ) 



(1 - 2a + a 2 ) d 2 + (ta - 2ia 2 + ia 3 ) d 3 + Q(d£) 



(3.55) 



(3.56) 

(3.57) 
(3.58) 
(3.59) 

(3.60) 



It is interesting that this sytematic expansion gives yet another first order correction to the non linear polarization 
- but of course it only differs from Brabec and Krausz's 0] and Porras's [I[ by terms like (1 — a) d T , which are second 
order corrections. Note again that Brabec and Krausz work in the case a = 1, whereas Porras allows a ^ 1, but does 
not consider nonlinear processes. 



IV. APPROXIMATIONS: SEWA, SEEA, AND GFEA 



The full generalised few-cycle equation (|3.48p has a rather complicated prefactor for the polarization term, and also 
the "extra" Trhs term. If we want to make approximations that reduce it to Brabec and Krausz's SEWA (slowly 
evolving wave approximation), Porras'sQ SEEA (slowly evolving envelope approximation), or some other form then 
we need to consider these terms in detail, and understand in what limits these terms might be simplified or neglected. 
Note that setting T~r — is a matter of chosen convention, and is not an approximation of any kind. 

For example, Brabec and KrauszQ have a discussion in their PRL on p3284, after their eqn.(8) about the various 
criteria that need to hold for their SEWA to be valid, all of which are that various quantities must be slowly varying 
as the pulse propagates along £. It is also instructive to see what parameter values or what terms are neglected to 
see how my generalised equation (|3.48|) reduces to Brabec and Krausz's SEWA (slowly evolving wave approxima- 
tion), Porras [4|'s SEEA (slowly evolving envelope approximation), or to my GFEA (generalised few-cycle envelope 
approximation), which is in some sense equivalent to a "best of" combination of the two others. 

1. SEWA: Brabec and Krausz[3j eqn.(6). In eqn. (13.48j) . set a = 1, use Tx = (since uj € = (x>o)and ignore Trhs'- 



qd ( A(>r ± ,£,T) 



2lTT 

Vi^(rj_,C, r) + -V (1 + id T ) S(rj_,e, r; AX-4.1) 



2. SEEA: PorrasQ eqn. (2). In eqn. (pHg)) . set B(r±,£,r; A) = 0, use Tx = and ignore T RHS : 
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gd e A(fx,£,r) = + tf)') A{? ± , £ r) + ' VlA^, £, r). (4.2) 

V Po / ^Po U + wC T J 

3. GFEA: Keeps the acuracies of both Brabec and Krausz .'5 and PorrasQ, whilst avoiding the more complicated 
parts of the full generalised cqn. (|3.47|) . In eqn. (|3.48|) . set 1 — a <C 1 and use Tx = and ignore Tuns'- 



qdtA(f> ± ,Z,T) = A(r ± ,Z,r) + ^ ^+^ B(r^,r;A) + ^ \ V 2 ± A(r ± , £, rX-4.3) 

V Po / Uq {1 + uj0 t ) 2{l + iaa T j 

The advantage of the generalised eqn. (|3.48|) over these is that I could easily put a = 1 — 6 and do an expansion 
for small S and go beyond both Brabec and KrauszQ and Porras^. 

The main point of complication with the generalised equation is the Trhs term, as it is second order in d T and 
c>£, as well as having higher order terms in d T from the diffraction contribution D' . These same terms (or simplified 
versions) were also neglected by both Brabec and Krausz [3j and Porrasj4j. Since the equation would be considerably 
easier to solve if Trhs were negligible, I will examine it in carefully in order to see what justification or constraints 
are required to do so. The term is: 



Trhs 
1 + icrd T 



-9! 



00 



ib' 



(4.4) 



Clearly this Trhs needs to be small compared to the other terms in eqn. (|3.48p if it is to be neglected. If it happens 
that it is small, then the other, non-negligible, terms sum to close to d^A(f±,^,T). So we can self-consistently ignore 
Trhs if the following condition holds: 



1 + iad T )~ 



Po 



<\dsA(? ± ,£,T)\ 



(4.5) 



A nice way of dealing with the presence of the d T terms is to fourier transform into the freqency domain, where 



id T — > n. 



(4.6) 



This enables us to avoid speculation about the possible time derivatives of A, and instead constrain its frequency 
components. However, this assumes knowlege of the all-time behaviour of the terms under consideration, so when 
being careful it might be better to use a time-windowed transform or similar. This could then give us constraints valid 
over some finite timescale relevant to the dynamics, without having to deal directly with instantaneous derivatives. 

The condition (|4.5|) can be broken into two parts, which are 



%A{?x,£,t)\ « \(l + iad T )d 6 A(f^,^T)\, 

2 



CTo 
00 



iD'\ A(f x ,Z,T) 



<£\dtA(r x ,£,T)\ 



(4.7) 
(4.8) 



The second of these I assume holds as a further consequence of the "first order" dispersion condition (|4.17p 
[SeeRef:dO] below. 



1. Note on the use of moduli 



The constraints we are attempting to apply are that the RHS term(s) has negligible effect on the propagation 
compared to that of the LHS term(s); the specific mathematical expression of this comparison is up to us. The 
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situation is complicated by the fact that either side can be complex, and will likely have a different complex argument 
(i.e. phase). Clearly the largest number we can make with the (hopefully small) LHS is given by the modulus, so that 
will give us an appropriate value for the LHS. What to do with the RHS is less obvious, because we would instead 
(to be cautious) want to pick a smallest reasonable value; but (e.g.) picking the minimum value of either the real or 
imaginary part would miss the point: indeed if the term were real, the smallest imaginary part would be zero; thus 
leading the condition to always fail. This leaves us with little choice (as far as I can see) but to use the modulus 
again; which in any case this would be the typical physicist's approach. 

Perhaps the best justification for applying the modulus to both sides in the comparison follows from the fact that a 
rotation in the complex plane applied equally to both terms should not affect the outcome. Hence we can rotate both 
terms so that the (hopefully large) RHS term becomes real-valued; now, since only the LHS might be complex, taking 
its modulus gives a useful upper bound on the significance of its contribution to the dynamics. This is equivalent to 
just taking the moduli of both terms. 



A. Evolution: 9| approximation 

I can now constrain the evolution of the pulse in £ by evaluating how to ensure that the <9| terms is negligible. 
Starting with eqn. (|4.7|) . 



\&A(? ± ,Z,t)\ « (l + ™)%i(fx,£,Q) 



(l + <7Lo)A(r ± ,Z,n) 
i(r±,£,0) 



(4.9) 
(4.10) 

(4.11) 

(4.12) 



because (i) I assume I can cancel (hopefully with no side effects) a derivative term from either side [SeeRefidO]; 
and (ii) I tighten the constraint somewhat by relying on er ~ 1 (not a ~ 1) and that only considering positive 
frequencies means that il > 0. Eqn (|4.7p is therefore revealed as a condition that the envelope function only changes 
slightly when propagated over distances £ ~ 1 (i.e. z ~ l//?o), and as such will depend on other terms in the evolution 
equation. 

More carefully, if we assume we know the £ behaviour of the envelope A, can use it to see how these constraints 
might hold in spatial- frequency (n) space (since £ <-> —in), 

\dlA(r ± ,£,T)\ « |(l + zad r )^A(r ± , K ,T)| (4.13) 

l 2 K 2 A(r±,K,T) 



cancelling i and k gives 
using moduli 



iK,A(r±, k, t) 
KA(r*±, k, r) 



ik (1 + io~d T ) A(r_j_, k, t) 
- (1 + iad T ) A(fj_,K, t) 
(1 + md r ) A(r*±, k,t) 



(4.14) 
(4.15) 
(4.16) 



This gives us an expression rather like that as if the propagation gave us different wavelength (wavevector) to that 
specified by the carrier exponential. We want the envelope to propagate with its spatially (£) behaviour to be peaked 
around small values of k, so that the bulk of the spatial variation is included in the carrier exponential." 

Note that the constraint will always be violated somewhere if the spatial bandwidth of the propagation is too large 
(even if it has an e.g. exponential fall off). The approximation therefore amounts to ignoring such violations of 
the constraint in the (spatial-frequency) wings of the propagation, on the basis they are "negligible" (which indeed 
seems perfectly reasonable). In a simulation, the LHS and RHS can be calculated as it progresses, and the frequency 
counted, and possible significance assessed (I have done this in OPA simulations). 

SeeNote:dO: There may be complications I do not see, because these terms (£>', etc) do contain derivative terms. 
However, since SVEA treatments ignore such issues, I do also. 



The spatial evolution sub-condition (|4.12p . the dispersion sub-condition (|4.17|) (c.f. condition (14. 8JI ); and hence the 
total condition (|4.5|) will only hold if the other non-negligible terms in eqn. (|3.48|l are similarly small compared to 
\A(r±,£,T)\, viz: 
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ViA(fx,£,r) 



2p$ (1 + zct<9 t ) 
2i7T (l + z<9 T ) 2 



(1 + iad T ) 



B(r ± ,Z,r;A) 



<C |A(f±,C,r)| 
< \A(r±,£,T)\ 

<C |A(f x ,e,r)| 



(4.17) 
(4.18) 

(4.19) 



Note that we could use the alternative eqn (|3.50p in condition (|4. 19|) . 

I will assume condition (|4.17|) implies that (|4.8p also holds [SeeRef:dO], leaving us with four conditions in total 
(|47T2l [47T71 [47T81 [47T9|) . Note that Brabec and KrauszQ] split ([4~T7]) into multiple pieces, which will be discussed in 
the following subsection. 



B. Dispersion: d T Approximation 



I will now treat the dispersion condition (|4.17p in the above approximations and associated conditions. Note that 
Brabec and Krausz[3] claim that their SEWA "does not explicitly impose a limitation on the pulse width"; however 
this is rather misleading as shortly afterward they introduce a pulse duration r p which is used in the inequalities 
constraining the material parameters, which then give the region of validity of the SEWA. Even the weakest statement 
we might make about t p needs to state that it does constrain the SEWA, because it (further) constrains the material 
parameters! For a given set of material parameters, there will be some pulse width limitation, although it might well 
be the few-cycles we hope to describe. 

I now break up eqn. (|4.17p into parts containing single factors of 7 m . Note that I need to exclude any terms 
including (3q and /3i in j m as they were treated separately in the analysis. I thus write 7^ = 7 m for m > 2, and 
7o = <xq, 71 = a.\. The inequality is then 



ui m <y' 



(3 m 



« \A(r ± ,tr)\ 
A(r±^,^) 



(4.20) 
(4.21) 



where the second line has been fourier transformed in time. Note that while both loq and 7' have units, does not, 
as it is the counterpart of the dimensionless (scaled) r (see (|4.6p ). 

To Do: Now get from to a t$ or Brabec and KrauszQ r p , by some physical motivation justifying t$ — 2tt 

I could treat the time-domain condition (|4.20p qualitatively as follows: introduce Tg ~ A/ \d T A\, the time for which 
a rate of change of d T A would accumulate (in absolute value) to A; and is thus some measure of how much time it 
takes the envelope to change significantly, and is thus hopefully something we can relate to the pulse width (or at 
least the width of one "bump" on the pulse envelope). Also, we assume that the higher derivatives of A are related, 
with \d™A/ml \ /A ~ r^ m [SeeRef:xl]. This is essentially the same parameter as Brabec and Krausz[3f's t p , but 
scaled into my dimensionless picture: t p = ts/luq. Using this t$ we can rewrite eqn. (|4.20p as 



I 1m 



< 1 



'5 - 



> O 1 ~ L~> 



(4.22) 
(4.23) 



This final condition is the same as those for a m and (3 m in Brabec and Krausz[3j, since (/?o and wavelength comment), 
and (L 7jm justification cf BK). 

SeeNote:xl:: The factoring of the ml into <9™ A seem the most mathematically most sensible thing to do by analogy 
to the terms in Taylor expansions, and expansions of exponential functions (etc). 

Note: Is it possible to derive a true time-domain treatment of the constraints by replacing the 7 m by <9 w fc(w), and 
relate it back to the time domain e(i — to) (or moments thereof). 

Note: the eqn. D' 2 terms. 
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C. Diffraction: approximations 



Treating the d T in condition (|4.18[) by fourier transform, as above, we can use the fact that for gaussian beams with 
a beam waist wq, we have \7 2 ± A(r±, £, r) ~ Wq 2 ; similar statements could be made for other typical beam profiles. 
This leads to the diffraction constraint on the SEWA becoming 



(1 + oft) /Sgtfljj > 1. (4.24) 

Comparing this to the comparable condition in Brabec and Krausz|3j after their eqn (8) (i.e. in my units PqWq 3> 1), 
we see that they are the same except for the fi term, so that my condition is in fact somewhat less restrictive than 
theirs. This is because considering positive frequencies only means that f2 > 0; and f2 ~ 1 for d T modulations much 
less than the variation of the carrier frequency. 



D. Nonlinearity: B(r±,£,T; A) approximations 



The nonlinearity constraint (|4.19|) is very complicated, so I ignore the term in square brackets [...]. Then, treating 
the d T in condition (|4.19|) by fourier transform, as above, and using a ~ 1, the nonlinearity constraint on the SEWA 
becomes 



2tt (1 + ail) 



n% (1 



> 



B(f ±) t,n;A) 



A(f±,£,n) 



(4.25) 



Comparing this to the comparable condition in Brabec and Krausz[l] after their eqn (BK8), we see that the same 
comments as for the diffraction hold - my condition is somewhat less restrictive than theirs. Of course it may be 
convenient to simplify the LHS of eqn. (|4.25p with various small Q expansions. 

In both cases I could instead include the d T corrections qualitatively by replacing d T — > t^ 1 before proceeding from 
(|4.18p and (|4.19j) : but even with t$ ~ 1 this will not alter the contraints greatly. 
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